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A self-contained approach to studying the unitary evolution of coupled qubits is introduced, 
capable of addressing a variety of physical systems described by exchange Hamiltonians containing 
Rabi terms. The method automatically determines both the Weyl chamber steering trajectory and 
the accompanying local rotations. Particular attention is paid to the case of anisotropic exchange 
with tracking controls, which is solved analytically. It is shown that, if computational subspace is 
well isolated, any exchange interaction can always generate high-fidelity, single-step controlled-NOT 
(CNOT) logic, provided that both qubits can be individually manipulated. The results are then 
\ applied to superconducting qubit architectures, for which several CNOT gate implementations are 

^ . identified. The paper concludes with consideration of two CNOT gate designs having high efficiency 

^v^j ' and operating with no significant leakage to higher-lying non-computational states. 
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I. INTRODUCTION 



Controllability of quantum mechanical systems has been the subject of numerous investigations in the last several 
r—v years Q, 0, 0, 0, 0, 0, 0, 0, 0, [IS Gil] • An important contribution by Khaneja et al. on time-optimal control 0,0] has 
led to the development of rf-pulse sequences for NMR-spectroscopy with nearly ideal performance [HI]. In Refs. [|[ 
II ' and it was assumed that the local terms in the Hamiltonian can be made arbitrarily large, which would allow an 
almost instantaneous execution of singlc-qubit operations. However, such hard control mechanism is not applicable to 
quantum computing architectures based on superconducting Josephson devices, in which the relevant computational 
^ ' subspace must be kept well isolated at all times. 

In this regard, the work of Zhang et al. on geometric theory of nonlocal two-qubit operations [HI, [HI] acquires special 
significance. The authors introduced a convenient, geometrically transparent description of SU(4) local equivalence 
\ classes and then used it to develop several implementations of quantum logic gates that did not involve hard-pulse 
control sequences. The description of entangling operations presented in Refs. [l3| and [l4[ is based on the fact H, 0] 
■^j- ' that any two-qubit quantum gate U € U(4) can always be written as a product, called the Cartan decomposition, 
00 ' 

|> ■ U = e^fcit/entfe, fci,A; 2 eSU(2)<g>SU(2), (1) 



with 

[/ cnt = e -(*/2H c l {r l ff 2+ c »<'l 2+ c 3 CT l '2,) - ( 2 ) 



in 
o 



The triplet of numbers c = (ci, 02,03) in Eq. ^ may be taken to represent the local class of U. In general, such 
representation is not unique due to the presence of symmetries mapping class vectors to other class vectors of the 
same equivalence class. However, it was shown in [l3| that the correspondence can be made unique if c is restricted 
to a tetrahedral region of 5R 3 , called a Weyl chamber. One such chamber is chosen to be canonical. It is described by 
the following three conditions [HI, [IB] : 

(i) 7T > Cl > c 2 > c 3 > 0, 

(ii) ci + c 2 < 7T, 

(iii) if C3 = 0, then c\ < tt/2. 

When a physical system evolves under the action of its Hamiltonian, c traces a trajectory inside the Weyl chamber, 
which explicitly shows the (continuous) sequence of dynamically generated local equivalence classes. For example, 
the Hamiltonian H = gafa^/^ generates the straight line c{t) = (gt, 0, 0) in the Weyl chamber. In this case, Eq. (fTJ) 
reduces to 

U{t) = e~ itH ee U ent (t), fci,fc 2 = l. (3) 
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After steering for a time £cnot = tt/2<? the system hits the gate 



U(t 



CNOT ) 
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71 



/ 1 -i\ 

1 -i 

-i 1 

V-i 1 



with 



7T/2X (1,0,0), 



(4) 



(5) 



belonging to the controlled-NOT equivalence class. By flanking £7(icNOT) with additional local rotations K\ and K 2 , 
any gate in that class can be made. For example, to make the canonical CNOT gate we can take 



CNOT 



e lTr/ V 



i(7T/4) a \ _i(7T/4)K- 



(6) 



Ki 



K 2 



When Rabi terms are present in the Hamiltonian, the steering trajectory is no longer a straight line. In Ref. 
the trajectory c(t) was calculated using the relation between the class vectors and the local invariants [l6j]. That 
method was applied in Ref. [ijj to a CNOT gate design for flux-qubits with SQUID-based controllable coupling. 

In the present paper we propose an alternative approach to finding the steering trajectory that does not rely on 
local invariants. Our goal is to develop a systematic procedure for calculating the entangling part U eni (t) of the 
time-dependent gate U(t) together with the accompanying it local rotations fci(t) and k 2 (t), so that the Cartan 
decomposition JT]) could be determined at every step of system's evolution. It turns out that due to a special property 
of the relevant to our problem generators of su(4) the closure under commutation and the existence of a central 
element — the local rotations required to implement JT]) can be chosen in a particularly simple form, which mimics 
the form of the local Rabi parts of system's Hamiltonian. Due to such simplifying form of k\ and k 2 , the full problem 
of steering can be analytically solved in the experimentally important case of tracking control. 

In the mathematical portions of this paper we will use the notation that is convenient for Lie algebraic manipulations: 

X k = (z/2K, XX = (i/2Kof , YY = {i/2)a\al, ZZ = (i/2)afa|, YZ = (i/2)**a$, ZY = (i/2)of<7», (7) 

with k — 1,2. The operators listed above form the Lie algebra Lq = span {Xx , X2 , XX, YY, ZZ, YZ, ZY} C su(4) 
whose commutators are given in the following table: 





X! 


x 2 


XX 


YY 


ZZ 


YZ 


ZY 













-ZY 


YZ 


-ZZ 


YY 


x 2 











-YZ 


ZY 


YY 


-ZZ 


XX 























YY 


ZY 


YZ 











-x 2 


-Xx 


ZZ 


-YZ 


-ZY 











Xx 


x 2 


YZ 


ZZ 


-YY 





x 2 


-Xx 








ZY 


-YY 


ZZ 





Xx 


-x 2 









(8) 



Later on, in sections devoted to applications, we will revert to the usual notation. 

Notice that it is possible to generate Lie algebras isomorphic to © by replacing the local operators (Xx, X 2 ) with 
either (Yx,Y 2 ) or (Zx, Z 2 ), without any change in our results. An example of this will be given in Sec. IIII C 31 



II. GENERAL CONSIDERATIONS 



Let us consider a generic time-dependent Hamiltonian 

iH(t) = n lx (t)Xx + n 2x (t)X 2 + g xx (t)XX + g vy (t)YY + g zz (t)ZZ + g yz (t)YZ + g zy (t)ZY, 



(9) 



whose scalar functions will be called the steering controls, or control parameters. The solution to the Schrodinger 
equation 



dU(t) 
dt 



-iH(t)U(t), 17(0) = 1, 



(10) 



3 



is a time-dependent operator U(t) € exp(Lo) C SU(4), which can always be written in the form [T| 



U(t) = e - Q ( t )^i-/ 3 ( t )^2 e -ci(t)XX-ci(t)YY-ca(t)ZZ e -C(t)Xi_-((t)X a 

fel(t) l/ent(«) M*) 



(11) 



The functions appearing in the exponents of Eq. (|11[) will be collectively referred to as the steering parameters, while 
the triplet (cx(t), 02(f), c^{i)) will be called the class vector, as usual [1 51 ] - In what follows, the class vector will be 
allowed to evolve on the full Cartan subalgebra Ac = span {XX, YY, ZZ} C Lq C su(4) rather than within the Weyl 
chamber, since projecting it onto the Weyl chamber can always be easily performed [151 ]. It is important to remember 
that at any given time t the choice of the steering parameters is not unique. Therefore, additional requirements (such 
as smoothness, initial conditions, etc.) must be imposed on the corresponding functions in order to determine the 
experimentally meaningful trajectory. 

Differentiating (fTTj) with respect to the time t gives: 



dU(t) 
dt 



= -[a'Xi + (3'X 2 + c[XX + e~ aXl e- pX2 (c' 2 YY + c' 3 ZZ) 



JX 2c aX! 



_ e - a X le ~pX 2e ~c 2 YY e -c 3 ZZ (( . /Xi + £ X2) e 



e 

. ZZ o c 2 YY o 0X 2o aX x] y/ f x 



(12) 



Here, each of the nested similarity transformations represents a rotation by some angle in a certain two-dimensional 
subspace of the Lie algebra Lq. For instance, 



3 c 3 zz 

e -C2YY YZe c 2 YY 



e 6 X\e 



-aXi 



YZe 



aXx 



X\ cos C3 + Y Z sin C3 , 
Y Z cos C2 + X2 sin C2 , 
YZ cos a + ZZ sin a, 



(13) 



etc. Using (fTUj) to perform algebraic manipulations in (fT2"j) and equating the resulting coefficients of the corresponding- 
generators on the right hand sides of (fTU|) and (fT2"|) , we get a nonlinear system of seven first-order differential equations 



1 
10 
10 
Ai A 2 
A 2 A x 
A 3 -Aa 
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-A 3 C 3 + A4C4 
A 4 C 3 - A3C4 
A1C3 + A 2 C 4 





C-2 
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9yz 








. 9zy . 



(14) 



M 

where the new variables 

C\ = cos c 2 cos C3 , C 2 = sin c 2 sin C3, C 3 = cos c 2 sin C3, C 4 = sin c 2 cos C3 , 

and 

Ai = cos a cos/3, A 2 = sin a sin/3, A3 = cos a sin (3, A 4 = sin a cos /3, 
have been introduced. Notice that 

det M = cos 2 c 2 — cos 2 C3. 

For simplicity, we choose 

ci(0) = a(0) = /3(0) = c 2 (0) = c 3 (Q) = C(0) = £(0) = 

to satisfy the initial condition U(0) = 1. 

The first equation in (I14p integrates immediately, 

ci(i) = I dr g xx (r), 



(15) 
(16) 
(17) 
(18) 



(19) 
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while the remaining system can be inverted to give 



a' 
& 

C 



^2x 



gvv(A 3 C33+AiC22)-gzz(A3C22+AiC33)-gyz(AlC33-MC22)-gzv(AlC22-A2C33) 

gyy(A3C22+A 4 C33)-gzz(A3C33+A i C22f-g I yz(A 1 C22-A2C3 3 )-gzy(A 1 C33-A2C22) 

dct M 

g yy Ai + g zz A 2 + g yz A 3 + g zy A± 
g yy A 2 + g zz Ai - g yz A A - g zy A 3 

-gyy(A 3 C3+AiCi)+gzz(A3C4,+AiC3)+gyz(A 1 C3-A2Ci)+gzy(.A 1 Ci-A2C3) 

-gyy(A3C i +A i C3)+gzz(A3C3+A i 8fj+gyz(A 1 C i -A2C3)+gzy(A 1 C3-A2C i ) 

dctM 



(20) 



where 



C 22 = cos c 2 sin c 2 , C 33 = cos c 3 sin c 3 . 
To make further progress, wc impose some restrictions on the form of the steering Hamiltonian. 



(21) 



III. ANISOTROPIC EXCHANGE WITH TRACKING CONTROLS 

Tracking (l9l | refers to steering with control parameters having the same enveloping profile defined by some function 
7(t). Notice that any time-independent Hamiltonian describes tracking with ^(t) — 1. Here we are interested in 
Hamiltonians 

iH(t) = j(t)piXi + n 2 X2 + gi(t)XX + g 2 YY + g 3 ZZ], (22) 

where gi(t) is a function of time and fix, f2 2 , g 2 , <? 3 are some constants. [It is possible to choose g\{t) arbitrarily because 
XX is central in L .] 



A. Solving the tracking control case 



Under these conditions, 



ci(t) 



dr 7(r)ffi (r) 



The remaining steering parameters will be found using the ansatz 

a(t) = C(t), /?(t)=£(t), 

or, equivalcntly, 

U(t) 



e -a(t)X 1 -(3(t)X 2 e -c 1 (t)XX-c 2 (t)YY-C3(t)ZZ e -a(t)X 1 -0(t)X 2 



(23) 



(24) 



(25) 



This ansatz works only for Hamiltonians given in (|22[) . For more general systems, another trick or numerical simula- 
tions based on (fTl"j) and (f2"0"|) should be used. 
The resulting system is 



a 



a' 
(3' 



= 7 



Qi 



(,g2A4-g3A 3 )C22 + (g2A3-g 3 A4)C33 

dct M 

q , (32^3-93^4)022 + (32^4-33^3)033 
1 i2 + dct M 

32^1 + 33^2 

g 2 A 2 + g 3 A 1 

-{g2A 4 ,-g3A3)C i -(g 2 A3-g3A l )C3 
dct M 

-(32^4-33^3)03 -(32^3-33^4)04 
dct M 



(26) 



The four equations for a' and (3' give 

7 [f2i(l + cosc 2 cosc 3 ) — J7 2 sinc 2 sinc 3 ] 



a 



(cosc 2 + cosc 3 ) 2 



7 [f2 2 (l + cosc 2 cosc 3 ) — Oi sinc 2 sinc 3 ] 
(cosc 2 + cosc 3 ) 2 



(27) 
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which determine a(t) and (3(t) after 02(f) and 03(f) had been found. Also, 

, _ (^1.93 - ^292) sin c 2 - (O x g 2 - ^253) sin c 3 (f2ig 3 - ft 2 # 2 ) sinc 3 - (Qi# 2 - ^233) shic 2 

A i — 9V7 ; s > A 4 — 



(9% -g|)(cosc 2 + C0SC3) 
The equations for c 2 and c 3 give 

(4±^ 2 -~ 2 

Using 



(9|-3l)(cosc 2 + cosc 3 ) 
(c 2 ±c 3 r = 1 \g 2 ±gz) 2 (A 1 ±A 2 f. 
(A X ±A 2 ) 2 = 1-(A 3T A 4 ) 2 , 



we get 



(c' 2 ±c 3 ) 2 = 7 2 



/ , \2 ( fi iT^) (sinc 2 ± sinc 3 ) 
(g 2 ±g 3 ) - - 



(cosc 2 + COSC 3 ) 2 

After applying the sum-to-product identities and re-arranging the terms, we arrive at 

d „:„ f^(t)±0 3 (t)\\ 2 , /7W./Z ^„N2 , .n^„ ; „^ 2 (f)±C 3 (f) 



By making substitution 



V 2 



( ff2 ± fl 3) 2 + (fiiTn 2 rsiii 



7(i) 



(52 ± 53 : 



/±(t) := sin 



c 2 (t)±c 3 (f) 



we can solve the resulting equation 



df±(t)Y , / 7 (*) 



^V(3 2± 53) 2 + ("iT^2)7±WJ = ^(52 ±53) 



by inspection. It is easy to see that 



f±(t) 



ff2 ± 33 



(g 2 ±g 3 ) 2 + (n lT tt 2 ) 2 



: sin 



( 32 ±53) 2 + («lT^2f 



j(t)cIt 



solves (f34| subject to (|18p. which together with (|23|) . (|25|) . (|27|) solves the tracking control case: 

C2,3(t) = arcsin (/+(*)) ± arcsin (/-(*)) ■ 



(28) 

(29) 
(30) 

(31) 

(32) 
(33) 
(34) 

(35) 
(36) 



B. Controlling the flow on the Weyl chamber 

Let us now assume that g\ is tunable, but otherwise independent of time. Then, given an experimentally realizable 
tracking mechanism j(t), a point on the Weyl chamber (alternatively, in the full Cartan subalgebra Ac) whose XX 
coordinate is c\ can be reached after steering for a time t\ satisfying 



7(r)dr = — . 

gi 



(37) 



This reachability condition is necessary, but not sufficient. Since the point is specified by a class vector c = (c\, c 2 , c 3 ), 
we have yet to determine whether the remaining coordinates c 2i3 can be realized by adjusting the Rabi frequencies 
rii 2- Here we will restrict our attention only to the points belonging to the XX axis, and thus having c 2j3 = 0. It is 
easy to see that in this case we must have 



(ci,0,0) 

1.2 



1 



lixngx 

Cl 



(52 - ff3) 2 ± 



27rmgi 



Cl 



(92 + 93) 2 



(38) 
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provided the integers n, m are chosen in such a way as to make the Rabi frequencies real. 

We can now write down a general condition under which a coupled qubit system directly generates controlled-NOT 
class corresponding to c — tt/2 x (1,0,0): 



n 

2g~i 



7 (r)dr = — , tt™ OT 



1 



v/(47i 9l ) 2 - (<7 2 - g 3 f ± y/(4m gi ) 2 - (<? 2 + <? 3 ) 2 



(39) 



Other approaches to CNOT gate design have been considered in Refs. 20, 21, 22, 2 



C. Tracking control of Josephson phase qubits 

1. Capacitive coupling with rf bias of Cliaf type 

In the rotating wave approximation (RWA) [24[ the dynamics of two resonant capacitively coupled phase qubits 
53, US HE H3, El, HI is described by the Hamiltonian 



Hi(t) = fr(t)/2)[Sl l0 f + g(afa% + o\o»)], g > 



(40) 



The Rabi term represents the action of an rf bias current applied to one of the qubits. It turns out that keeping just 
one such local term suffices to generate controlled-NOT logic [HI]. The condition = 53 = gives f+(t) = f-(t), 
which leads to 



c\{t) = g \ 7(r)dr, ci(t) = 2 arcsin 
Jo 



and 



ait) =£li I dr „ 

/ I + COSC2 



, /3(t)=0. 



, c 3 {t) = 0, (41) 



(42) 



The time-dependent gate is therefore 



U(t) = e-( i /2)aTr e -(V2)(ci<T^J+c 2 <T«<T») e ~(i/2)aa^ 



which becomes an element of controlled-NOT class, provided [IE 



/■*CNOT 

/ 1 (r)dr = n/2g, fi™ OT = 9y/W - I 
Jo 



(43) 



(44) 



with n = 1, 2, 3, 

We may use Result 1 of Ref. [13] to state the following applicability condition for the RWA: 
The solution to the Schrodinger equation with the RWA Hamiltonian (|4"U|) approximates the solution with exact H 
(reduced to computational subspace; see Ref. [l8[ for details) in the sense that if Qi/uj 1 (weak perturbation) and 
lu = e (resonant condition), then || - 0rwa(O — V ; oxact(0ll = 0(VLi/uS) whenever < t < 0(ui/Qi). Here, u) is the bias 
frequency, e is the computational level splitting. For UCSB architectures [3l| with qubit coupling g ~ 10 GHz and 
level splitting u> < 100 MHz, ftf NOT /w - 10~ 2 . 

For calculations that go beyond the RWA in the context of Josephson phase qubits coupled to nanomechanical 
resonators see Ref. [33 . 



2. Inductive coupling with rf bias of Qiaf + tyP e 



For inductively coupled qubits [Tfl, |33|, |3J, |35|, |36|, |37|, |38|, |39j driven by local rf magnetic fluxes the Hamiltonian in 
the RWA is [Hi 



H 2 (t) = (7(t)/2)[n x a? + n 2 o% + g{pfoi + a\a\ + ka{a z 2 % g > 0. 



(45) 
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Using (|36p . the steering trajectory is found to be 



ci(*) = 9 / 7(t)c!t, 



C2,3(i) = arcsin 



1 + fc 



± arcsin 



1 - fc 



: sin 



(|V(i + fc) 2 + Pi-^)/g] 2 ^ 7 (r)d 

sin f - fc) 2 + [(fii + ^ 2 )/5] 2 / 7M<* 

V(l-fc) 2 + [(«i + ^2)/5] 2 \2 VV ' /yJ io 



where 



= e -( i /2)(a^r+/3<T2) e -( 4 /2)(ci^ C rJ+ C2( T« CT «+ C 3 -j ! ( Tj) e -(i/2)(a ( Tj : +/3^)^ 

with a and /3 calculated from Eq. (|27[) . The CNOT class is generated by setting [19l ] 

f^CNOT 



v/(4n) 2 - (1 - fc) 2 ± v/(4m) 2 - (1 + fc) ; 



(46) 



(47) 



(48) 



For example, for g — 1.00, fc = 0.10 and n = m = 1, the Rabi frequencies are fii = 3.8716, H2 — 0.0258. The 
corresponding Weyl chamber steering trajectory for j(t) — 1, with parameters measured in units of it/2, is shown in 
Figures [111 and [31 



3. Inductive coupling with dc bias of Q,i{o~l — <r|) type 

Because of the {X\,X^) — > (Yi,!^) — > {Z\,Z-i) "symmetry" mentioned in Sec. HI it is possible to devise an 
alternative CNOT implementation based on the Hamiltonian for inductively coupled qubits acted upon by dc fluxes: 



H 3 (t) = (7(*)/2)[Oi(o? - a 2 2 ) + g{ka\a z 2 + of of + afaf ) 



(49) 



The effect of such bias is to "move" system's energy levels by equal amounts in opposite directions (the process known 
as detuning). One important feature of this implementation is that for any \k\ < 1/2 it is always possible to generate 
controlled-NOT logic by choosing Rabi frequencies < |f2i|/<? < 1. This is important when perturbation is required 
to be small (see Section [V] for a more general approach). 
We have, 



ci(t) = kg 7(r)dr, c?.{t) — c${t) = arcsin 
Jo 



V 



1 sin ( g^/l + ^/g) 2 / j(t)cIt] 

l + (fii/<?) 2 V Jo J 



and 



a(t) = -f3(t) = Y^ dT 



COS 2 C2 



(50) 



(51) 



where the steering parameters (a, (3, c\, C2, C3) are now associated with the operators (Z%, Z%, ZZ, XX, YY). The 
time-dependent gate is given by 



U(t) — e _ ^/ 2 ^ Q ( CT l _CT 2) e -(V2)(ciCTj ! CT2+C2(o-iCT2>CT^CT*)) e -(i/2)a((Tj ; -tTj)^ 



which implements CNOT class, provided 

P^CNOT 



7 (r)dr = n/(2kg), fif NOT = g\J (2kn) 2 — 1, 



where (2fcn) 2 > 1. Several examples of this implementation are listed in Table HI 

Figures |U [5] show the steering trajectory for g = 1.00, fc = 0.10, and the Rabi frequency f2i = 0.6633. 



(52) 



(53) 
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TABLE I: Generation of controlled-NOT logic with c = 


(•7t/2) x (1,0,0) using inductively coupled flux qubits 


subject to 


symmetric dc detuning (fii/2)(crf — erf). The Hamiltonian 


is given in (149[) . Here, -y(t) — 1. 




k icNOT, units of 7r/2gr 


n 


«1 /<? 


0.100 10 


6 


0.6633 




7 


0.9798 




8 


1.2490 


0.050 20 


11 


0.4583 




12 


0.6633 




13 


0.8307 




14 


0.9798 




15 


1.1180 


0.025 40 


21 


0.3202 




22 


0.4583 




23 


0.5679 




24 


0.6633 




25 


0.7500 




26 


0.8307 




27 


0.9069 




28 


0.9798 




29 


1.0500 



IV. DISCUSSION 



We now discuss limitations and possible extensions of the proposed method. 

The most significant limitation comes from restricting the local terms to form a homogeneous pair (such as, for 
example, (Xi,X 2 )). By adopting such restriction we were able to isolate a special subalgebra L of su(4), given in 
Eq. ([8]), that contains a central clement. The fifteen-dimensional problem was then reduced to a nonlinear system 
of "only" seven first-order differential equations, one of which completely separated from the others. By making a 
certain ansatz, the analytical solution in the tracking control case has been found. 

We can extend this approach to Hamiltonians with arbitrary combinations of Rabi terms, such as (X%, Y2), etc. 
The dimensionality of the problem would increase, but it would still be possible to write down and solve — most 
likely, numerically — the corresponding system of differential equations. 

For Hamiltonians containing homogeneous local terms with arbitrary time dependence the following useful ansatz 
can be identified: 

Case 1. For 



H(t) = (l/2)[Q 1 (t)af+g 1 (t)o-^+g 2 (t)afa : 



2J- 



use 



63(f) =0, /?(t)=£(t)=0, 
which corresponds to the Cartan decomposition 

Ult) = e -( ,; / 2 )«^ e -(V2)(ci^<Tf+c2<^) e -(V2)C^ 

Eq. (PO)) then reduces to 

S7i — (?2 sin a cos C2 / sin C2 
<?2 cos a 
<?2 sin a I sin C2 



'a 1 ' 








c 





(54) 
(55) 
(56) 

(57) 



Case 2. Anisotropic exchange with symmetric detuning. This case generalizes the detuning Hamiltonian considered 
in Section [III C 31 by allowing arbitrary time dependent controls: 



H(t) = (l/2)[ni(t)(of - <rf) +9i(t)afo-* 2 + fl2 (tK<7f + g 3 (t)c-fo-%) 



(58) 



In this case we use 



a (t) = -/3(<), at) = -m, 



(59) 
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U(t) = e _ (V 2 )"( CT l- CT !) e -(V 2 )( c l CT l cr 2+ c 2C r rc r 2 E + c 3 cr l '2) e -(V 2 )C(o-l-0"2) > 



Eq. (|20|) now becomes 



\a'- 

















^i — (.92 + 33) cos a sin a(cos C3 sin C3 — sin c 2 cos C2)/ (cos 2 ci — cos 2 C3) 

(72 cos 2 a — 33 sin 2 a 
(73 cos 2 a — 32 sin 2 a 
(32 +33) cos a sin a/ (cos C2 sin C3 + sin C2 cos C3) 



Case 5. For systems described by 

H(t) = (1/2)PM*K + na(*)of +fli(*Kaf + 52 (t)(«7M + af^)] 



c 2 (t)=c 3 (t), C(t)=0, 



E7(t) = e -( i /2)("CT 1 ' ; +/3^) e -( l /2)[ Cl( T^ CT J+c 2 ( CT » CT «+ CT J ( T|)] e -(i/2)C^ 



corresponding to 



Notice that in this case we cannot use Eq. (I20|) directly because matrix M is not invertible, as can be seen 
(JTTj) . Instead, the original system (fT4]) has to be re- written in accordance with the constraints imposed by 
then get 



(60) 

(61) 

(62) 
(63) 

(64) 

from Eq. 
(1631 . We 



1 























1 







c 2 




& 




n 2 








A x + A 2 


-(A 3 


- A 4 )C 3 




c' 2 




32 








A 3 -A 4 




f A 2 )C 3 . 




-C - 




. 



with the variables defined as before, and 



The system can now be inverted to give 



\a'- 








c' 2 









Mi 



det Mi — cos c 2 sin c 2 . 



^1 + 32 (A3 — A 4 ) cosc 2 / sinc 2 
2 + 32 (A 3 - A 4 ) sinc 2 / cosc 2 
32 (Ai + A 2 ) 
-32 (A 3 - A4)/ cos c 2 sin c 2 



(65) 
(66) 

(67) 



which can be solved numerically. 



V. REDUCING LEAKAGE TO NON-COMPUTATIONAL STATES. 



Here we describe two controlled-NOT gate implementations satisfying certain constraints that must be imp osed on 
Josephson phase qubits [3l[ in order to make leakage to higher-lying (non-computational) states small [40l . l4ll . [42] | , 
while maintaining the high efficiency of the gate. The relevant conditions are: 

1. Hamiltonian: 

Hit) = (7(t)/2)[nxx<7? + fi»20f + n »W + n v^2 + + + 5«^2 + ^2 + (68) 



2. Coupling constants: 



3 > 0, \k\ < 0.5. 



(69) 
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3. Number of H applications: 



4. Rabi frequencies: 



N = 1. (70) 
\Cl xi \,\n yi \,\n zi \ <g, i=l : 2. (71) 



5. Efficiency: 



V ■= > 2.5. (72) 

9 %ate 

All these constraints can be satisfied by directly steering toward the target belonging to the CNOT equivalence 
class with entangling part C/ C nt(icNOT) represented by the class vector c = (tt/2) x (1,0,0). The canonical CNOT 
gate can then be made out of J7(fcNOT) = &it4nt(icNOT)&2 by performing additional local rotations K\ and as 
usual. 

The two implementations are: 

1. Symmetric dc detuning: In this case the Hamiltonian is 

<ldcW = (7(*)/2)[OiK-af) + n 2 ( f 7f-<7l)+ foK-qf) +ff«of + o?of + fcafol)], (73) 

symmetric dc detuning 

or, alternatively, 

<i. dcW = (7(i)/2)PiK + <%) + fiaK + <tf) + naK-al) +<?(afaf + ^ + fcafaf)], (74) 

symmetric dc detuning 

with 

^i=5, (75) 

where 7 (i) represents experimentally available tracking control, and the superscript w refers to the corresponding 
choice of the x and y Rabi parts. The relevant control parameters have been found numerically and are listed in Table 

m 

Figures[H]and[7]show the steering trajectory for g = 1.00, k = 0.050, with Rabi frequencies = 0.0133, ^3 = 0.7575. 

2. Asymmetric dc detuning: In this case the Hamiltonian is 

H { J m . dc (t) = ( 7 (t)/2)[fi 1 (of - o») + Q 2 (a$ - a\) + p 3 cr* - O^rf ) +g(af^ + a\a\ + kafa z 2 )}, (76) 

asymmetric dc detuning 

or, alternatively, 

dc(*) = (7(*)/2)[ni(of + af) + n 2 (a* 2 + o*) + [Sl 3 af - iltf) +g(afa x 2 + afa* + kafa^)}, (77) 

y v ' 

asymmetric dc detuning 

with 

n 1= n 3 = g. (78) 
The corresponding steering controls are listed in Table IIII1 



VI. CONCLUSION 



In summary, we have proposed a self-contained approach to steering on the Weyl chamber and applied it to the case 
of anisotropic exchange with tracking controls, which was solved analytically. It was shown that if architecture allows 
for local manipulation of individual qubits, any exchange interaction can generate CNOT quantum logic. The results 
were then used to identify several CNOT gate implementations for superconducting Josephson qubits, including the 
ones that are capable of suppressing leakage to noncomputational states without significant reduction in the gate's 
efficiency. 
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TABLE II: Generation of controlled-NOT logic with c = (tt/2) x (1, 0, 0) using inductively coupled flux qubits driven by weak 
local perturbations and subject to symmetric dc detuning (fi3/2)(crf — erf). The Hamiltonian is given in Eq. (|73[l or (174[l . 
where £li/g = 1. Here, 7(t) = 1. 



k 


tcNOT, units of 7r/2(? 


^ NOT /s 


^ NOT /s 


V 


0.000 


1.595776 


0.000000 


0.755502 


2.5066 


0.001 


1.595775 


0.000264 


0.755503 


2.5066 


0.002 


1.595774 


0.000529 


0.755505 


2.5066 


0.003 


1.595772 


0.000793 


0.755509 


2.5066 


0.004 


1.595769 


0.001057 


0.755515 


2.5066 


0.005 


1.595765 


0.001322 


0.755522 


2.5066 


0.010 


1.595731 


0.002644 


0.755582 


2.5067 


0.025 


1.595496 


0.006614 


0.756001 


2.5071 


0.050 


1.594657 


0.013257 


0.757500 


2.5084 


0.075 


1.593263 


0.019961 


0.760001 


2.5106 


0.100 


1.591321 


0.026758 


0.763506 


2.5136 


0.150 


1.585843 


0.040779 


0.773549 


2.5223 


0.250 


1.569080 


0.071908 


0.806036 


2.5493 


0.350 


1.547002 


0.111865 


0.856120 


2.5856 


0.450 


1.530753 


0.178169 


0.927506 


2.6131 


0.490 


1.550430 


0.240369 


0.966790 


2.5799 


0.493 


1.561200 


0.254105 


0.971189 


2.5621 



TABLE III: Generation of controlled-NOT logic with c = 


(tt/2) x (1,0,0) using 


inductively coupled flux qubits driven 


by weak 


local perturbations 


and subject to asymmetric dc detuning (Q.z/2)al — (f24/2)<rf 


. The Hamiltonian is given in Eq. (|76l) 


or J77JI, 


where £li/g = ^3/1 


/ = 1. Here, j(t) = 1. 








k 


tcNOT, units of n/2g 


«2 NOT /ff 


Of OT /9 


V 


0.000 


1.553771 


0.000000 


0.402539 


2.5744 


0.001 


1.553770 


0.000179 


0.402541 


2.5744 


0.002 


1.553768 


0.000358 


0.402548 


2.5744 


0.003 


1.553766 


0.000537 


0.402558 


2.5744 


0.004 


1.553762 


0.000715 


0.402574 


2.5744 


0.005 


1.553757 


0.000894 


0.402593 


2.5744 


0.010 


1.553716 


0.001789 


0.402757 


2.5745 


0.025 


1.553430 


0.004475 


0.403902 


2.5749 


0.050 


1.552414 


0.008974 


0.407988 


2.5766 


0.075 


1.550736 


0.013523 


0.414781 


2.5794 


0.100 


1.548418 


0.018150 


0.424259 


2.5833 


0.150 


1.541995 


0.027780 


0.451143 


2.5940 


0.250 


1.523410 


0.050016 


0.535559 


2.6256 


0.350 


1.501442 


0.081649 


0.659439 


2.6641 


0.450 


1.488962 


0.141937 


0.826279 


2.6864 


0.500 


1.515587 


0.220268 


0.938373 


2.6392 


0.506 


1.539498 


0.251771 


0.959755 


2.5982 
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FIG. 1: Steering trajectory generating CNOT class in the case of rf-biased inductively coupled flux qubits, Eq. (|45[) . 
g = 1.00, k — 0.10, fii = 3.8716, Q.2 = 0.0258, 7(t) = 1. The steering parameters are given in units of 7i"/2. 



Here, 




FIG. 2: Local rotation a accompanying the steering trajectory shown in Fig. [T] 




FIG. 3: Local rotation /3 accompanying the steering trajectory shown in Fig. [T] 




FIG. 4: Steering trajectory generating CNOT class in the case of inductively coupled flux qubits subject to dc symmetric 
detuning, Eq. f|49[) . Here, g = 1.00, k — 0.10, Sli = 0.6633, y(t) — 1. The steering parameters are given in units of tt/2. 




FIG. 5: Local rotations accompanying the steering trajectory shown in Fig. [4] 
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FIG. 6: Weyl chamber steering trajectory generating CNOT class in the case of inductively coupled flux qubits driven by weak 
local perturbations and subject to symmetric dc detuning, Eq. (|73[) . Here, g = 1.00, k = 0.050, O2 = 0.0133, f^3 = 0.7575, 
-y(t) = 1. The steering parameters are given in units of tv/2. 



